QUANTUM FLUCTUATIONS OF CHARGE AND PHASE TRANSITIONS OF A 
LARGE COULOMB-BLOCKADED QUANTUM DOT 
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We analyze ground-state properties of a large gated quantum dot coupled via a quantum point 
contact to a reservoir of one-dimensional spinless electrons whose interactions are characterized by 
the Luttinger liquid parameter g. We find that the classical step-like dependence of the average 
number of electrons on the dot n as a function of the gate voltage n x is preserved under certain 
conditions, the presence of quantum fluctuations notwithstanding. We point out that in its low- 
energy limit the problem is dual to that of a single-junction SQUID with Caldeira-Leggett dissipative 
coupling, and analogous to the classical problem of multilayer adsorption, and so is related to the 
classical statistical-mechanical problem of a one-dimensional Ising model with exchange interactions 
decaying as the inverse-square of distance. This Ising universality class further subdivides into what 
we term (i) the Kondo/Ising class and (ii) the tricritical class, (i) For systems of the Kondo/Ising 
class, the n(n x ) dependence is always continuous for g > 1, while for g < 1 (repulsive electrons in the 
constriction) the n(n x ) dependence is continuous for sufficiently open dots, while taking the form of 
a modified staircase for dots sufficiently isolated from the reservoir. At the phase transition between 
the two regimes the magnitude of the dot population jump is only determined by the properties 
of the reservoir and given by p 1//2 . (ii) For systems in the tricritical class we find in addition an 
intermediate regime where the dot population jumps from near integer value to a region of stable 
population centered about a half-integer value. In particular, this tricritical behaviour (together 
with the two dependencies already seen in the Kondo/Ising class) is realized for non-interacting 
electrons, g = 1. 

PACS numbers: 71.10.Pm, 73.43, Jn, 73.21.Hb, 72.10.Fk 



I. INTRODUCTION 

Quantum dots are nanometer-scale structures hosting 
a few to a few thousand electrons and coming in a 
variety of forms. Quantum dots can be realized as small 
metallic islands Q , gated two-dimensional electron gases 
in semiconductor heterostructures [0-pl, and most re- 
cently as short segments of metallic carbon nanotubes || . 
The latter two realizations are of particular interest as 
they exhibit a remarkable tunability. Through the appli- 
cation of metallic gates upon either the heterostructures 
or upon the nanotubes, the chemical potential together 
with the strength of any putative dot-lead coupling can 
be tuned. The associated flexibility afforded to the ex- 
perimentalist has helped to make quantum dots and their 
related phenomena an active area of research. 

One such phenomenon that has been closely studied, 
both experimentally and theoretically, is Coulomb block- 
ade. Coulomb blockade in its simplest form is a fascinat- 
ing effect arising from the discrete nature of electrons 
combined with their mutual repulsion j?| . The origin of 
the effect is entirely classical. Imagine charging a dot of 
capacitance C biased by a gate voltage Vq- The classical 
electrostatic energy of the dot, E c i = Q 2 / (2C) — QVq, is 
minimized if the dot charge Q equals CVq. But the dot 
charge can only change in the increments of the electron 
charge e, and the dimensionless gate voltage n x = CVq/ e. 
is generally not an integer Jqj. Typically there is only 
one integer closest to n x which determines the number 
of electrons on the dot to be n = Q/e = [nj. Half- 



integer values of n x play a special role because now two 
states of the dot with n = n x ± 1/2 electrons have the 
same energy. Thus as n x is varied between two near- 
est half-integers, the dot charge is fixed and quantized 
(the Coulomb blockade) but at half-integer values of n x 
the blockade is momentarily lifted, and the dot charge 
changes discontinuously by e, leading to the celebrated 
Coulomb staircase 0. 

This simple classical picture is modified by quantum 
charge fluctuations arising from the coupling to a reser- 
voir of electrons. Quantum mechanical tunneling of 
charge between connecting leads and the dot has been 
argued to smear the Coulomb staircase even at zero tem- 
perature j^]. Although originally treated perturbatively 
H , the smearing finds its strongest support in a connec- 
tion of the quantum dot problem to a highly anisotropic 
Kondo problem Jicj ]. In the limit of low transparency, i.e. 
a small tunneling matrix element, and of large Coulomb 
repulsion, the charge dynamics on the dot can be ap- 
proximated by that of a two-level system, with the lev- 
els corresponding to charge states of the dot differing by 
unit charge p0| . Changing the charge on the dot then 
becomes akin to flipping an impurity spin- 1/2. In Kondo 
language the net charge on the dot parallels the magne- 
tization of the impurity spin while the deviation of the 
gate voltage from the half-integer n x is akin to a Zeeman 
field applied to the Kondo spin. For a single-channel 
Kondo problem, the spin susceptibility is well-known to 
be finite at zero applied field. This correspondingly ne- 
cessitates the smearing of any discrete charge jump with 



1 



the variation of the gate voltage. 

The connection to a Kondo problem applies in the 
weak-tunneling limit. The opposite high-transparency 
limit was treated by both Flensberg and Matveev 
PH - (For a comprehensive review of both cases, see Q; 
for more recent work on this limit see 0].) On the ba- 
sis of this work it was possible to conjecture that the 
smeared Coulomb staircase of the weak-tunneling limit 
smoothly evolves as tunneling is increased into a strictly 
linear function of the gate voltage. These predictions 
p2| found partial confirmation in recent experimental 
work Jl5| ]. The agreement between theory and experi- 
ment was limited by finite temperature effects, with the 
temperature of the quantum dot being far in excess of 
the putative Kondo temperature. 

In this paper the role of quantum fluctuations is re- 
examined for a large quantum dot - the dot whose charg- 
ing energy e 2 /(2C) is substantially larger than the dis- 
tance between the energy levels on the dot. We show 
that as the strength of the tunneling is tuned, changing 
the dot from open to closed, a modified Coulomb stair- 
case reappears under certain circumstance with tunnel- 
ing non-zero. In particular we find that quantum fluc- 
tuations need not destroy the Coulomb blockade. We 
in part arrive at our general conclusion by relating the 
zero-temperature dynamics of spinless electrons to the 
classical statistical mechanics of a one-dimensional Ising 
model with ferromagnetic interactions decaying as the 
inverse-square of distance. Within the Ising universal- 
ity class, we find that all experimental systems must be 
either of the Kondo subclass (to which the connection 
to the Kondo system given in JTcj ] applies) or of a novel 
tricritical type. 

The organization of the paper is as follows. In Sec- 
tion II we set up the problem in general terms, showing 
how one understands the problem as belonging to the 
Ising universality class. We demonstrate how the two 
subclasses we identify can be differentiated in terms of 
the form of their effective potentials. In Sections III and 
IV we examine these two subclasses, the Kondo and tri- 
critical subclass, in detail. In Section V we discuss our 
results. 



II. FORMULATION OF THE PROBLEM 

The experimental setup we have in mind is a gated 
quantum dot placed in a large magnetic field (so that the 
electrons can be considered spinless) coupled to a single 
reservoir (lead) via a narrow constriction - the quantum 
point contact |12|] . A related setup has been experimen- 
tally realized . 

Regardless of the nature of the reservoir, the prob- 
lem is effectively one-dimensional ]ll],[l2|. If the reser- 
voir is high-dimensional, the point nature of the contact 
allows one to dimensionally reduce the electrons in the 



lead through a partial wave expansion, keeping only the 
s- waves, to a one-dimensional problem in much the same 
way the problem of Kondo impurities embedded in bulk 
metals can be reduced to a one-dimensional problem [ fl6| . 

The most general phenomenological Euclidean action 
for the problem is 



where 



E (if,,n x ) = U{<f>) 



2C 



(1) 



(2) 



Here r is the imaginary time, and 4>(t) is a fluctuating 
electron number field - the expectation value of </>(t) gives 
the average number of electrons on the dot, n —< <j) >. 

The first kinetic term in the above action arises in sev- 
eral steps. In the limit of large dot size the time scale 
by which electrons scattered through the point contact 
return to the scattering point is much larger than any 
other time scale in the problem 11 T^. The point con- 
tact then connects in effect two reservoirs of electrons. 
Upon dimensionally reducing the electrons in the higher- 
dimensional leads, one obtains a kinetic term appropriate 
to one-dimensional electrons. Bosonization of the elec- 
trons leads to a kinetic term of the form 



dxdrid 11 ^^). 



While the term in the action describing tunneling be- 
tween the dot and the lead is naturally restricted to a 
spatial point [first term of (|^)], it was shown jl2| that the 
Coulomb charging [second term of (|^)] can also be rep- 
resented as acting at a single spatial point. As these two 
interaction terms are restricted solely to a single point we 
can integrate out the bosonic degrees of freedom away 
from the point contact. Doing so reproduces the first 
term of (Q) with g = 1 which corresponds to noninter- 
acting electrons in the constriction |l7| . The frequency 
cutoff, A, that the kinetic term is equipped with is of or- 
der the ratio of the sound velocity to the electron density 
(both in the vicinity of point contact). 

In (Q) we look at a somewhat more general prob- 
lem by allowing the one-dimensional electrons to interact 
through the inclusion of the dimensionless Luttinger liq- 
uid parameter g which is generally different from unity. 
This case can be realized if both the reservoir and the dot 
are one-dimensional systems (having common g) such as 
quantum wires, metallic carbon nanotubes or fractional 
quantum Hall edges Jl8| . 

The perturbing energy Eo((j),n x ), as previously indi- 
cated, consists of two terms acting at a single point: a 
term by which electrons tunnel through the point con- 
tact and a term expressing the charging energy of the 
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dot. The former may generally be represented as a term 
both periodic and even in 0, 

U{4>+i) = u{4>); u{4>) = u{-4>). 

The latter, charging, term takes the natural form ^{4>~ 
n x ) 2 . The function U(<p) having its minima at integer and 
maxima at half-integer (f> reflects the preferred tendency 
for the dot to host an integer number of electrons. The 
classical limit Hj is recovered by "closing" the dot, i.e. by 
increasing the amplitude of U to infinity. On the other 
hand, the regime of perfect transmission is reached by 
"opening" the dot, i.e. by decreasing the amplitude of U 
to zero. While in the limit of almost perfect transmission, 
the form 

U ((/>) ~ cos 2-7T0 

was derived in jjJI through bosonization, the precise 
functional form of U ((/>) should depend on the details of 
the constriction. 

The above action finds applications beyond that of 
quantum dots. The inspection of (Q) and (§) reveals that 
in the low-energy limit the Coulomb blockade problem 
is dual to that of a single-junction SQUID |l9| with a 
Caldeira-Leggett |2(J dissipative coupling to the environ- 
ment. Here the magnetic flux through the SQUID loop 
is dual to the charge of the dot Q (the flux quantum $o 
is dual to the electron charge e), and so all our conclu- 
sions will have their dual counterparts. This example of 
electric-magnetic duality is discussed in further detail by 
Likharev @. 

Alternatively by viewing the imaginary time coordi- 
nate r in (|lj) as a fictitious space direction we arrive at 
an effective Hamiltonian defining a classical statistical 
mechanics problem pi] . The above action can be then 
recognized as describing a one-dimensional version of the 
problem of multilayer adsorption p2|-^|. In this prob- 
lem a film in contact with a liquid is adsorbed onto an 
attractive substrate. The film thickness (controlled by 
external pressure) is finite as the free energy of the ma- 
terial of the film exceeds that of the liquid. For thick 
films the competition between the substrate attraction 
and the confining effect of external pressure can be mod- 
eled by a parabola centered near optimal film thickness 
p2| . For a crystalline film there is also an additional pe- 
riodic contribution into the free energy which oscillates 
with the number of adsorbed atomic layers. If the field 
<f> is regarded as a number of adsorbed atomic layers, 
then (^) can be identified with the uniform part of the 
film free energy. The \uj\ dependence of the first term of 
(|l|) translates into a 1/r 2 repulsion occuring between two 
parts of the film having different heights and separated 
by distance r. The dependence of the film thickness upon 
pressure is analogous to n(n x ). 

In what follows we will be interested in the n(n x ) 
dependence. This quantity has been shown to be di- 
rectly accessible experimentally. Through fabricating 



two quantum dots in close proximity, the authors of |15| ] 
were able to measure the charge on one dot by monitor- 
ing the current through the other. The n(n x ) dependence 
has an inversion symmetry about every half-integer point 
n = n x , and the whole dependence can be reconstructed 
from any segment of n(n x ) of unit length. Thus without 
the loss of generality we restrict ourselves to the interval 
of n and n x between and 1. 
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FIG. 1. Schematic plots of the ground-state energy as a 
function of the dot occupancy n and corresponding n(n x ) de- 
pendencies for case (i) (Ising/Kondo) as discussed in the text. 
The first three rows in the figure represent successively larger 
values of n x beginning with n x — 1/2 in the topmost row. 
Equilibrium dependence of the dot population n on the gate 
voltage Tlx is presented by solid curves in the sketches of the 
bottom row. The hysteretic parts of the n(n x ) dependencies 
are shown by dotted lines, and the arrows on the hysteresis 
loops reflect the direction of the change of n x . 

In order to compute the n(n x ) dependence we need 
to minimize the full ground-state energy E(n,n x ) of the 
system, in principle computable from (Q) and @ . Quali- 
tatively the outcome can be anticipated without calcula- 
tion by noticing that akin to Eo((f>,n x ), the ground-state 
energy E(n, n x ) as a function of n must be a sum of a pe- 
riodic and parabolic functions. Here there are two main 
cases to consider. 

Case (i): In the first case the maxima of the periodic 
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part of E are sufficiently sharply peaked and the ground- 
state energy E as function of n in the interval < n < 1 
sees either one or two minima. Fig.l illustrates this case. 
The top row shows two possible sketches of the ground- 
state energy as a function of n at the degeneracy point 
n x = 1/2, and two subsequent rows demonstrate how 
these curves tilt as n x increases. The solid curves at the 
bottom row show equilibrium n(n x ) dependencies found 
by minimizing the ground-state energy. 
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FIG. 2. Schematic plots of the ground-state energy and cor- 
responding n(n x ) dependencies in case (ii) (tricritical) . The 
plot is organized as in Figure 1. 

In the first scenario, pictured in Fig. la, there is a single 
minimum in E, the position of which smoothly changes 
as n x is varied. The population of the dot, n(n x ), is 
then continuous and will look as sketched in the bottom- 
most figure of Fig. la. In the presence of two minima, 
pictured in Fig. lb, a modified Coulomb staircase is re- 
covered upon adiabatic variation of n x . This staircase 
differs from its classical counterpart in that the pop- 
ulation plateaus will have a nonzero slope and the jump 
in dot population occuring at n x = 1/2 between the two 
degenerate minima is less than unity. The inter-plateau 
jumps are first-order phase transitions, and as such they 
can be accompanied by hysteretic phenomena. Indeed 
if one changes n x non-adiabatically, the system stays at 
the left well until the well ceases to exist and only then 
will it jump into the remaining right well. The result- 



ing n(n x ) dependence will deviate from the equilibrium 
staircase. Upon decreasing n x a different path, n(n x ), 
will be taken and a hysteresis loop will be produced as 
shown in the bottommost sketch of Fig. lb by the dotted 
lines with arrows. 

These two scenarios are distinguished primarily by the 
amplitude of the periodic part, and one can imagine going 
between the staircase and continuous n(n x ) dependencies 
by closing or opening the dot. Alternatively however, one 
can imagine varying the reservoir interactions (varying g) 
or changing the dot size, i.e. its capacitance C . 

Case (ii): If the maxima of the periodic part of E are 
sufficiently flat, the ground-state energy, E, as a func- 
tion of n can have as many as three distinct minima. 
Fig. 2, in a fashion similar to Fig.l, illustrates this more 
complicated case. The three cases presented in Fig. 2 are 
connected by varying the amplitude of the periodic part 
of the energy. 

In the column (a) of Fig. 2, the amplitude of the pe- 
riodic part is weak and the ground-state energy has a 
single minimum for any n x . This results in a continuous 
dependence n(n x ) depicted at the bottom of Fig. 2a. As 
the amplitude of U is increased, one arrives at column 
(b) of Fig. 2. When the gate voltage is tuned to de- 
generacy point, n x = 1/2, the ground-state energy has a 
global minimum at n = 1/2 and degenerate metastable 
minima at n close to and 1. If one increases n x adia- 
batically away from 1/2, the number of electrons on the 
dot changes continuously (from a half-integer value) until 
the depths of the "half-integer" and the right "integer" 
minima become equal. This happens at a critical value 
of n x in the range, 1/2 < n x < 1, and a first-order phase 
transition onto the "integer" plateau occurs. Upon fur- 
ther increasing n x , the state of the dot changes continu- 
ously until the next "half-integer" minimum comes into 
play. As a result, a staircase with both integer and half- 
integer plateaus is produced, as sketched at the bottom 
of column (b). We will arrive at qualitatively the same 
n(n x ) dependence if "integer" minima are not present at 
the top energy curve but develop and become global in 
tuning the gate voltage away from n x — 1/2. 

As the amplitude of the periodic part of the potential 
is increased further one arrives at column (c) of Fig. 2. 
Here the ground-state energy has degenerate minima at 
n close to and 1 and a metastable minimum at n = 1/2. 
If initially the dot state corresponds to the leftmost well, 
and one starts adiabatically increasing n x away from 1/2 
the energy curve tilts, and the system jumps into the 
rightmost well. Upon further increase of n x the average 
number of electrons on the dot increases continuously 
until the next degeneracy point is reached. As a result 
a modified Coulomb staircase, as given at the bottom of 
column (c) in Fig. 2, is recovered and the presence of the 
third metastable minimum does not qualitatively affect 
the dot occupancy (compare with Fig. lb). Hysteresis is 
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also possible in this case as shown in Figs. 2b and 2c. 

We can reconsider both case (i) and (ii) from a unified 
viewpoint by focusing solely on the charge degeneracy 
points n x = 1/2. Qualitatively Eq(4>, 1/2) appears as in 
the top rows of Figs.l and 2 and can be described by 



E 



-u 2 s 



U 6 S 



with s = 0—1/2 (the coefficients 1*2,4 can be of both signs 
while uq > 0). When employed in (|l]), we arrive at a 
one-dimensional continuous-spin Ising system with 1/t 2 
ferromagnetic exchange interactions. A small deviation 
from the charge degeneracy points can be modeled by 
applying a weak fictitious magnetic field h ~ n x — 1/2, 
i.e. by adding a —hs term to the expansion (an equivalent 
argument using the language of the Kondo model has 
been given by Matveev fl(i| , |l2|| ) . If U4 > 0, the sixth order 
term can be ignored and we obtain the standard long- 
ranged Ising model, corresponding to case (i) discussed 
above. With U4 < stability demands we keep uq, and 
we obtain a long-ranged version of the tricritical Ising 
model pf| and so reproduce the findings of case (ii). 

We note that recently, while studying the problem of 
occupation of a resonant level (for example, an impurity 
or a small quantum dot) coupled to a Luttinger liquid, 
Furusaki and Matveev |^6| also noticed a connection to 
the classical one-dimensional Ising model with interac- 
tions decaying as inverse square of distance. 

In the next two sections we consider the cases (i) and 
(ii) in more quantitative detail. 



III. KONDO/ISING TYPE 

In this section we consider a dot whose potential, 
Eq (0) , can be approximated by a two- well potential and 
is thus in the ordinary (as opposed to tricritical) long- 
ranged Ising universality class. The truncation of charge 
fluctuations of the quantum dot to that of a two-level 
system performed in Jl0| maps the dot problem to the 
conventional Ising model with long-ranged 1/t 2 interac- 
tions pTj] or equivalently the single-channel Kondo model 

The physical consequences of this equivalence can be 
derived from an examination of the renormalization- 
group flows of the long-range Ising/Kondo model. To 
understand these flows we follow the treatment of the 
dual dissipative SQUID problem p{|, and introduce two 
parameters, the distance between the two wells of the po- 
tential q, and the (dimensionless) inter-well rate of tun- 
neling A. Under a renormalization-group transformation 
these parameters flow at lowest order as follows @,|9| : 



dA 



drn(Ai) 



(1 - ?7s)A, 



dq 2 



dln(At) 



-g 2 A 2 . (3) 



Here t is the running scale and the equations describe 
how the two parameters renormalize upon integrating out 



high-frequency modes. The initial conditions, Ao and 
qo, can be computed from the form of the action in (Q). 
An example of such a computation done in the SQUID 
context may be found in p9| . 

The flow diagram corresponding to (||) is sketched in 
Fig. 3. 
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FIG. 3. The flow diagram corresponding to (|3|). The sta- 
ble part of the A = fixed line is shown in bold, and the 
arrows indicate the direction of the flow. The locus of the 
free-electron models, g = 1, is presented schematically by the 
dotted line. 

The flows can be divided into two main regions, corre- 
sponding to the ordered and the disordered phases of the 
Ising model. 

In the former case are systems satisfying q 2 /g> 1 and 
A§/2 < q 2 /g - 1 - \n(q 2 /g), i.e. those to the right of 
the separatrix A 2 /2 = q 2 / g — 1 — ln(q 2 / '<?), are carried 
by the flow to a stable fixed line, A — 0, q 2 / g > 1. 
The absence of tunneling between the two minima of the 
potential indicates that we are in the ordered phase of the 
Ising model. The functional form of n(n x ) is as shown in 
Fig. lb, and the dot population discontinuity, An (given 
by q), at n x = 1/2 satisfies the inequality 



An> q 1 ' 2 . 



(4) 



The equality is reached at the Ising phase transition 
point, Aq/2 = q 2 /g — 1 — \n(q 2 /g), where the popula- 
tion jump An = g 1 / 2 depends only on the properties of 
the reservoir. This universal relationship is analogous to 
that between the magnetization jump, the amplitude of 
inverse-square interaction and the phase transition tem- 
perature of the Ising model |p^-|30[| . 

Although systems to the left of the separatrix A 2 /2 = 
q 2 1 g — 1 — ln(g 2 / g) are taken by the flow outside of the 
range of applicability of Eqs. (||), the effective tunnel- 
ing, A, between the two minima of the two- well potential 
can be seen to increase. We are thus in the disordered 
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phase of the Ising model. Here the dependence, n(n x ), 
is continuous as in Fig. la. Although the Ising transi- 
tion is continuous, in going from Fig. la to Fig. lb, the 
n(n x ) dependence changes discontinuously. In the disor- 
dered phase, the ground-state energy slightly away from 
n x = 1/2 is given by 

E = (7rh/g£) < s > 2 -h < s >, 

where the dependence on the correlation length £ is dic- 
tated by the one-dimensional nature of the problem. 
Minimization with respect to < s > leads us to the equa- 
tion, 



(a) = n - 1/2 ~ (gt/(2irh))(n x - 1/2), 



(5) 



describing the shape of the smeared population step. 

In classical limit of closed dot, Ao = 0, the distance 
between the minima of the two-well potential is unity, 
go = 1, which implies that if the electrons in the con- 
striction are repulsive, g < 1, these systems belong to 
the ordered phase of Ising model. Allowing weak tunnel- 
ing, i.e. sufficiently small Ao, renormalizes the popula- 
tion step downwards, and we arrive at the n(n x ) depen- 
dence of the kind shown in Fig. lb. At sufficiently large 
critical Ao c the Ising phase transition takes place, and 
for Ao > Ao c the dependence n(n x ) is continuous, as in 
Fig. la. 

In this same classical limit, Ao = 0, but with the elec- 
trons in the constriction attractive, g > 1, we find our- 
selves on the unstable part of the A = fixed line. Here 
allowing infinitesimally small tunneling, Ao, immediately 
destroys the classical staircase, thus leading to the con- 
tinuous n(n x ) dependence, Fig. la. For g > 1, the flow 
starts out vertically away from the A = line. In order 
to compute £ to lowest order in Ao, the renormalization 
of q given by the second of Eqs.(||) can be ignored. Inte- 
grating the first of Eqs.(||) up to a scale at which A « 1 
determines the correlation length to be £ w A" 1 Aq^ 1 a ' . 
Physically the difference between the repulsive and at- 
tractive cases can be understood by noticing that inter- 
particle attraction enhances quantum fluctuations in the 
constriction, thus making it easier for the electrons to en- 
ter and leave the dot; the effect of repulsion is opposite. 

The free-electron case, g = 1 is the most delicate to 
treat because it separates regimes which are stable and 
unstable against perturbations of infinitesimally weak 
tunneling. If the function U ((f>) is a periodic comb of 
delta-function wells of large but finite amplitude, then 
one still has go = 1- Solving Eqs.(^) with the initial con- 
dition qo = 1 and Ao -C 1, we find £ rs A -1 exp(7r/2Ao). 
If this £ is inserted in (JsJ) , then the shape of the smeared 
population step to leading order in Ao coincides with 
that given in Ref. jl(J. For all other periodic functions 
U(4>) one has go < 1) an d the locus of the corresponding 
models is shown in Fig. 3 by the dotted line. If the initial 
values go and Aq are not below the outgoing, leftward 



separatrix in Fig. 3, then £ w A -1 exp(ci/A ) where the 
constant c± interpolates between ir/2 (the q = 1 axis) and 
unity (the left separatrix) . If the initial values are below 
the separatrix, then one has £ w A^ 1 exp[(c2/Ao) In Aq -1 ] 
where C2 is another constant. To conclude, if the elec- 
trons in the constriction are noninteracting and we are 
in the Kondo/Ising universality subclass, then the n(n x ) 
dependence is always continuous, as shown in Fig. la, for 
arbitrarily small tunneling coupling . 



IV. TRICRITICAL ISING TYPE 

The results derived so far examine the case where at 
n x = 1/2 the effective potential energy can be approx- 
imated by a two-well potential. But this assumption is 
invalid for potentials whose periodic part has sufficiently 
flat maxima and thus whose energy dependencies appear 
qualitatively similar to the top row of Fig. 2. Moreover, 
if even the two-state assumption is correct for the bare 
potential, renormalization of truncated two-state poten- 
tial by zero-point motion might be qualitatively different 
from renormalization of the full potential. In such cases 
the tricritical behavior as pictured in Fig. 2 will be over- 
looked. 

Therefore it is important to reexamine the problem 
without making the two-state assumption. In this sec- 
tion we consider exactly such a case where approximat- 
ing the potential energy as two wells is inappropriate. In 
what follows we choose the periodic part of the potential 
entering (|^) to be U (4>) ~ cos 2iz(j), and rewrite (Q) as 

E ((f>,n x ) = —[irb{4> - n x f - ^-cos27r<£]. (6) 

Here b = e 2 g / (2-nC Ah) is the dimensionless charging en- 
ergy, while the dimensionless parameter v in the limit 
v <C 1 can be expressed in terms of the reflection ampli- 
tude at the point contact jl2| . This potential (^) models 
Coulomb blockade [p"2|, S QUIDs fl9[| , and multilayer ad- 
sorption phenomena |22]-|24||. Surprisingly we find that 
the model (g) does not belong to the Kondo/Ising uni- 
versality subclass although naively for n x — 1/2 the two- 
state assumption would appear to be reasonable. 

This model is readily understood if the classical limit, 
g — ► 0, is taken in the constriction. Then the fluctua- 
tions of the field, <ft, controlled by the first term of (|l|) 
are frozen. To find the n{n x ) dependence, we have to 
minimize (|J) with respect to <j> = n. The outcome is 
known from the SQUID context [^9|: for v > b the n(n x ) 
dependence qualitatively looks like that shown in Fig. lb 
(hysteretic SQUID regime). If, on the other hand, v < b, 
the n(n x ) dependence is of the kind presented in Fig. 
la. The phase transition at v — b is continuous with the 
population jump which grows from zero as one enters the 
staircase phase. 



6 



A. Scaling analysis 

To understand the role of quantum fluctuations in the 
constriction (nonzero g) we first treat Eo(tp,n x ) pertur- 
batively (i.e. v small, putting us in the high-transparency 
limit). The corresponding lowest-order renormalization- 
group equations 



vb 9 ~ 
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1 



dv 



d\n{M) 



= (l-<?)v, 



fib 



dln(At) 



= b, 



(7) 



parallel those given by Huse |23|] in the context of the 
adsorption phenomena. The parameters v and b of (|^) 
serve as initial values to Mj. We note that due to nonan- 
alytic | a; | dependence in (Q) there is no renormalization 
of g to any order. This is in contrast to the adsorption 
phenomena where the corresponding renormalization of 
an interface stiffness coefficient is an important part of 
physics. 

With b = 0, the first equation in (^) is known to de- 
scribe the physics of a Luttinger liquid in the presence of 
a weak point inhomogeneity p7| : there is a phase tran- 
sition at g = 1 which separates insulating, g < 1, from 
conducting regimes, g > 1. The second equation, solely 
due to the scaling transformation, tells us that the charg- 
ing energy is always a relevant perturbation. In g < 1 
regime, v does not decrease under the flow, and Coulomb 
blockade can potentially occur. In this regime integrat- 
ing the first equation of (0) until umI provides us with 
a length scale w A -1 !; 1 '^ 9-1 ^, i.e. a longitudinal cor- 
relation length |3l[] . In the absence of the charging term, 
at scales exceeding £„ , the insulating behavior sets in and 
the fluctuations of the field, </>, are suppressed. 

The presence of the charging energy introduces another 
length scale into the problem, w (A6) _1 , as can be 
inferred either from (fi ) and (g) , or from integrating the 
second equation of (0). In the absence of the periodic 
term at the scales exceeding £b, the fluctuations of <fi & r e 
suppressed. 

The evolution of the Coulomb blockade as a function 
of the system parameters can be understood qualitatively 
as a result of the interplay between the scales £„ and £b 
(both assumed to be much larger than the microscopic 
scale A -1 ). If £5 < the fluctuations of are sup- 
pressed by the charging term before the periodic term 
comes into play. This implies that the physics is domi- 
nated by the charging term, the dot population n tracks 
the gate voltage n x with electronic discreteness a small 
effect. The n{n x ) dependence is continuous as in Figs. la 
and 2a. If, on the other hand, > £„, the fluctuations 
are suppressed by the periodic term before the charging 
effects come into play. This means that the physics is 
dominated by the particle discreteness with n(n x ) corre- 
spondingly in the form of a staircase. The two regimes 
are separated by a phase transition occuring at ~ £„, 
i.e. for 



ln(l/v) 
ln(l/6)' 



(8) 



We note that (||) has roughly the correct limit b w v as 
g^O, and implies that for g = 1 a phase transition oc- 
curs at some v which is 6-indcpcndcnt in the limit b — ► 0. 
We now show that the scaling condition in (|^) under- 
spccifies the physics and that there are in fact two phase 
transitions occuring when ~ £„. 



B. Variational analysis 

For a more quantitative analysis we use Feynman's J3S 
variational principle for the ground-state energy: 



E < E x = E + (T/h) <A-A Q >o 



(9) 



where T is the temperature, and h/T has a meaning of 
the system size in the r direction; the T = limit will be 
taken at the end. The notation <>o denotes expectation 
values computed using an arbitrary reference action Aq>, 
and Eq is the ground-state energy corresponding to Aq. 

This method has been remarkably successful in analyz- 
ing the roughening phase transition p3| ], multilayer ad- 
sorption phenomena |22| , wetting transitions [ p4| , and a 
variety of problems of quantum mechanics and quantum- 
field theory |35j . Perhaps most relevant to our purposes, 
has been the application of the variation principle to 
the problem of quantum Brownian motion in a periodic 
potential ||(| which is mathematically identical to the 
single-impurity problem |T^| and formally the b — limit 
of the Coulomb blockade problem ([l]) , (0) . Here the vari- 
ational method provides exact description of the phase 
transition at g = 1 Mj. 

It is physically reasonable to select the trial action Aq 
in a Gaussian form similar to that in |B2l and |36l : 



Ao=—[J ^M|0M| 2 + Am J 'dr^-nf]. (10) 

-A 

The dimensionless variational parameters which include 
the familiar number of electrons on the dot n and a new 
parameter m (controlling the extent of fluctuations about 
n) are selected to minimize E\ in (|^). 

Using ( pij| ) the upper bound E\ entering @ can be 
computed. Ignoring unimportant additive and overall 
factors we have 

E x - (1 + b) ln(l + m) - bln(m) 

+ -1 [27r 2 &(n - n x ) 2 - v(—^—) 9 cos 27ml. (11) 

1 + m 

Minimizing E\ with respect to m and n we arrive at the 
system of two equations 



m — b = v( r cos 27m 



(12) 
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(13) 



1 —v cos 2nn ' 



(16) 



which are mathematically identical to those of Weeks ]22[| 
derived in the context of adsorption phenomena. 

The most valuable feature of the variational ap- 
proach is its nonperturbative nature. The perturbative 
renormalization-group treatment (|?|) is included as a spe- 
cial case as can be seen from an argument which parallels 
that of Brezin, Halperin, and Leibler originally given in 
the context of wetting transitions |57| : 

Let us split the fluctuating number field <f> into a 
"slow", 0o, and a "fast" part, (pi, so that 4> = 4>a+4>i, and 
"integrate out" fast degrees of freedom <j>\ treating (0) as 
a perturbation. This turns into its renormalized (but 
unrescaled) counterpart 

E L = —[Trb(cf> - n x f - V AI^L1 cos2w( f>], (14) 
g 2vr 

where the subscript is dropped for brevity, and at the 
scale L all high-frequency fluctuations are integrated out. 
We note that equations (0) and (|l^ ) are equivalent in 
terms of the physics they describe. The number of elec- 
trons n on the dot, and the correlation length L can be 
found from the conditions similar to those of |37||, 



[dE L /d<f>] n = 0, (2TrAh/g){LA)- 1 = [d 2 E L /d<f> 2 



(15) 



where the dependence on L in the second equation is 
dictated by the \u>\ dependence in (|]). It is straightfor- 
ward to verify that substituting (|lj) in (|lj) produces two 
equations which are just the m <C 1 limit of (|l^) and Jl3| ) 
provided the identification m — (LA) -1 is made. 

To summarize, the m«l limit of the variational anal- 
ysis is equivalent to the perturbative renormalization- 
group treatment. 

Below we analyze the coupled system of equations ( p"2] ) 
and (|l3|), and whenever multiple solutions are found, 
the one giving the lowest ground-state energy ( pi] ) is se- 
lected. We make two preliminary observations. Firstly 
we note that solving ( |l2| ) provides us with a periodic 
m(n + 1) = m{n) dependence which after substituting 
into (11) results in a function which is a sum of a peri- 
odic and parabolic functions as we already anticipated. 
Secondly, for n x integer or half- integer Eq.([l3j) has a so- 
lution n = n x as expected. 



C. Free electrons in the constriction: g = 1 



For (? = 1, Eq. (12) admits an analytic solution, and 
thus the ground-state energy (^) can be explicitly com- 
puted. The results are simplest in the limit m < 1. In 
this limit the solution to (O) is given by 



valid for b -C 1 and v < 1. Equation ( pif ) shows that 
fluctuations around integer dot population are always 
smaller compared to those for a half-integer number of 
electrons on the dot. Substituting (|l^) into ([l3]) we find 
the n x (n) dependence which parallels the result of Weeks 



v sin 2im 



2n(l — v cos 27m) 



(17) 



The remarkable feature of (|17| ) is its universality - there 
is no dependence on the dot capacitance in the I) « 1 
limit. We emphasize that the result ([l7]) goes beyond 
the original model (^j): any even periodic function U(4>) 
from (^J) can be decomposed into a Fourier series consist- 
ing of cosine functions. For g = 1 the leading harmonic 
[accounted for by (0)] responsible for ( |l7| ) is marginal 
[see the first of Eqs.([7M while higher harmonics are irrel- 
evant. In the limit b — > the contribution of higher-order 
harmonics vanishes (unless the corresponding Fourier co- 
efficients are unreasonably large which we assume to be 
not the case). 

At g = 1 the variational argument gives two criti- 
cal values of the dimensionless amplitude, v c i and v C 2'. 
the lower w cl marks the onset of the regime with half- 
integer plateaus while the upper w c2 gives the transition 
to a regime solely with integer plateaus, i.e. a modified 
Coulomb staircase. 

To determine v c \ one can study the monotonicity of 
the n x (n) function as given in (|l7|). For sufficiently small 
V, the n x (n) dependence is monotonic. For solving 
(|l7j ) to lowest order in v reproduces the result of Matveev 

n = n x — (v/2ir)sin27m x . 

As v becomes larger, n x (n) ceases to be monotonic. 
Defined when dn x /dn vanishes, this happens at v c \ = 
V3/2 < 1 and cos 2nn c = 1/V3 |||. The correspond- 
ing gate voltage n x is not a half-integer, a first indication 
that the lower critical value v c i separates continuous, Fig. 
2a, from discontinuous behaviour, Fig. 2b, with popula- 
tion plateaus centered both around the integer and half- 
integer number of electrons on the dot. 

To better understand the phase diagram we look at 
the ground-state energy. Combining (|ll| ) and ([l2]) for 
g = 1, with (|l6|), taking the limit m, b <C 1 and dropping 
unimportant additive constants we arrive at the form, 

Ei(n,n x ) ~ 6[ln(l - v cos 2im) + 2n 2 (n - n x ) 2 ]. (18) 

Comparing the periodic parts of Jl8| ) and (||) we con- 
clude that quantum fluctuations sharpen integer minima 
and flatten half-integer maxima. The minimization of 
Ei with respect to n reproduces (|l7|). For v < v c \ the 
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ground-state energy E\ always has a unique minimum 
which results in a continuous n(n x ). For v just slightly 
above v c \ — and n, x = 1/2 the ground-state energy 

still has a global minimum at half-integer dot popula- 
tion. But upon deviating from n x = 1/2, a near integer 
minimum develops and eventually becomes global. The 
phase transition at v c \ is continuous. 

As v is increased beyond v3/2, we come to the second 
transition. The approach to this transition is character- 
ized by the deepening of two metastable minima close to 
n = and n = 1 (Fig. 2b) with n x = 1/2. The criti- 
cal value v C 2 marks the point at which two near integer 
minima become equal in depth to the half-integer mini- 
mum. This is a first-order transition between the regimes 
sketched in Figs. 2b and 2c. The second critical value can 
be estimated accurately analytically from the condition 

Ei(l/2, 1/2) = ^(0,1/2), 

resulting in v c2 « tanh(7r 2 /4) « 0.9857 < 1. The true 
v C 2 is marginally smaller as the competing minimum is 
shifted slightly away from n — 0. For v C 2 < v < 1 
there are only population plateaus centered around inte- 
ger number of electrons on the dot, although the ground- 
state energy always has a meta-stable minimum corre- 
sponding to a half-integer dot population. 

Examples of the three regimes separated by v c x and 
v C 2 are plotted in Fig. 2 with values of v = 0.8,0.97, and 
0.988 for (a)-(c) respectively. The continuous segments 
of the n(n x ) plots were found by inverting ( [i"7| ) while the 
location of the first-order jumps was determined from the 
minimization of (|l8|). The hysteretic behavior in Figs. 2b 
and 2c arises when the inversion leads to multiple solu- 
tions and was determined by continuing the monotonic 
regime solution beyond the range where it is energetically 
favorable. 

It is possible to determine the slopes dn/dn x of the 
plateaus (given that that the slope is non-zero, the lan- 
guage of plateaus is used figuratively) . The slope of the 
integer plateaus given by 1 — v is rather flat in the range 
v c \ = v3/2 < v < 1. This flatness is a direct conse- 
quence of the sharpness of the near integer minima of 
the energy curves. For the half-integer plateaus over this 
same range of v we find the slope to be much larger, l+v. 
This steepness is a consequence of the shallowness of the 
half-integer minima in the energy curves. Alternatively 
one can recast this in the language of the tricritical Ising 
model pq ]. Two adjacent integer plateaus are then re- 
garded as two ordered Ising phases of opposite magneti- 
zation. The half-integer state of the dot is identified with 
the disordered (zero-magnetization) phase of the Ising 
model. On the n(n x ) curves of Fig. 2b the disordered 
phase coexists with both ordered phases. Equivalently 
upon increasing v, the disordered phase disappears grad- 
ually in the interval between v c \ and v C 2- On the other 
hand, for the systems of the Kondo/Ising type only or- 



dered phases can coexist on the n(n x ) curves. Corre- 
spondingly, increasing the amplitude of U in (0) beyond 
some critical value leads the disordered phase to disap- 
pear altogether. 

As was demonstrated in Section IIIB, the m«l limit 
of the variational treatment is equivalent to the perturba- 
tive renormalization-group analysis of the problem based 
on Eqs.(^) or equivalently Eqs.(|l4|) and (|l^). We specu- 
late on the grounds of the following plausibility argument 
that for the special case of free electrons, g = 1, such an 
analysis becomes exact in the limit 6^0: 

For 6 = the right-hand side of the first of Eqs.(^) is 
believed to vanish to all orders in v if and only if g = 1 
J36|. Equivalently this means that on taking the limit 
3 = 1 and b = in ( p~4| ) one would arrive at an exact 
renormalized potential. On the other hand, in deriving 
Eqs.©, ©, and © a small finite b was kept, and the 
6^0 limit was taken in the final formulas. It appears 
reasonable that higher order terms in b not included in 
( |l4| ) will not affect the results in the b — ► limit. Specif- 
ically, we argue that Eqs.©, @, and © solve the 
Coulomb blockade problem exactly for v < 1 and 6^0. 

For v > 1 the assumption m« 1 cannot be justified, 
and we need to return to the original system of equations 
(|il"|)-(|l3|). From numerical studies we find that for b <C 1 
and v > 1, the n(n x ) dependence rapidly evolves into a 
near perfect staircase. 

The foregoing analysis presumes b <C 1. As b increases, 
the slope of the integer plateaus, the critical values v c i 
and v C 2, and the range of existence of the intermediate 
phase with half-integer plateaus all become larger. How- 
ever no qualitative changes occur from what was found 
in the b — > limit. 

D. Interacting electrons in the constriction: g ^ 1 

In general we find that our conclusions at g = 1 are 
quantitatively but not qualitatively modified by the pres- 
ence of interactions. Interactions change the level of 
quantum fluctuations in the system, making it easier 
(g > 1) or more difficult (g < 1) for the electrons to 
enter or leave the dot. 

In the limit b,vV~ x < 1, Eqs. © and © give the 
continuous dependence 

n = n x — [vb 9 ~ l /2-k) sin 2i:n Xl (19) 

thus generalizing the g = 1 result of Matveev [jl2j . If the 
electrons inside the constriction are attractive, g > 1, 
Coulomb blockade effects are substantially weakened. 
The range of existence of the continuous n{n x ) depen- 
dence grows compared to the noninteracting case; the 
critical values v c \ and v C 2 grow similarly. As g goes to 
infinity the number of electrons on the dot n rapidly ap- 
proaches the gate voltage n x . If, on the other hand, the 
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electrons are repulsive, g < 1, the effects of Coulomb 
blockade are enhanced. The range of existence of the 
continuous n(n x ) dependence shrinks, and the critical 
values v c i and v C 2 decrease approaching their common 
v = b limit as g — > 0. 




FIG. 4. The phase diagram of the model (g) for b = 0.1 
and < g < 1. The phase with integer-half-integer sequence 
of the population plateaus exists between the lower, v c i(g), 
and upper, v C 2(g), critical curves. 

Fig. 4 shows the numerically derived phase diagram of 
the model (||) for b = 0.1 and < g < 1. In between the 
two curves for v c \ and v c i lies the phase with half-integer 
plateaus. These two curves meet at v c \ — v c i — b in the 
limit g — > as expected from (||). In further agreement 
with the scaling result (||) in the & < 1 limit, the critical 
curves take the form 

VdAg) = ai^b 1 ^ 9 

with a 1)2 (0) = 1, ar(l) = V3/2, and a 2 (l) w 0.98. 

We also verified that except for g — the n x = 1/2 
ground-state energy (|ll|) always has a minimum at half- 
integer dot population: g = 0, v = b is a tricritical point 
of the model (g). 

Fig. 5 shows the numerically derived n(n x ) dependen- 
cies for g = 1/2, b = 0.01, and (a) v = 0.13, (b) v = 0.16, 
and (c) v = 0.2. The values of v are selected to show con- 
tinuous (a), and discontinuous n(n x ) dependencies with 
integer-half-integer (b) and integer-integer (c) sequences 
of population plateaus. 
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FIG. 5. The dependence of electron number n on the dot 
on the gate voltage n x for g — 1/2, b — 0.01, and (a) v — 0.13, 
(b) v — 0.16, and (c) v — 0.2. Hysteresis loops are shown by 
the dotted lines and the direction of change of n x is indicated 
by the arrows. 
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V. DISCUSSION 

The most intriguing prediction of Eqs.(p])-(|l3|) is the 
existence of a range of parameters where the popula- 
tion of the dot jumps from a near integer value to a re- 
gion of stable population centered about a half-integer 
value. The origin of this effect is the persistence of the 
half-integer minimum in the population dependence of 
the ground-state energy at the charge degeneracy points. 
The effect is quantum - the bare energy (|6|) does not have 
this property |3S[| . For the effect to occur, the periodic 
part of the ground-state energy must have sufficiently flat 
maxima which the bare energy (JgJ) does not have. 

The effect is not entirely unexpected and resembles the 
classical effect of an inverted pendulum |39| : if the point 
of support of a plane pendulum in a uniform gravitational 
field is subject to rapid oscillations, a new minimum in 
the effective potential energy can be induced. The quan- 
tum fluctuations of the dot population resemble oscil- 
lations of the support, and the mechanical analog of the 
SQUID (and of the Coulomb blockade problem) is a slight 
modification of the simple pendulum 1 4^] . The contribu- 
tion to the pendulum's potential energy from gravity is 
given by the second term of (||) with 2ir(p analogous to 
the pendulum's angular position. The first term of (|^) 
describes a harmonic torsion bar one end of which is at- 
tached directly to the pendulum axle while the other is 
held at some fixed angle, this angle being analogous to 
the gate voltage, n x , of (||). The inverted pendulum anal- 
ogy is at best suggestive, and our effect is weaker - we 
only find a flattening of the maxima of the periodic part 
of the dot energy. 

Direct experimental observation of the intermediate 
phase with integer-half-integer sequences of plateaus 
might encounter several difficulties. The fluctuations 
around the half-integer dot population are large and 
likely a model- independent effect. Moreover in the 
regime of parameters, i « 1 and g < 1, the range of 
existence of the intermediate phase is extremely narrow. 

If the plateaus were to be observed, the corresponding 
experimental signal of entering this phase from the con- 
tinuum phase would be a divergence of the gate-reservoir 
capacitance (proportional to dn/dn x ) at the lower criti- 
cal value of the amplitude of the periodic potential, v c i, 
at two values of the gate voltage placed symmetrically 
about n x = 1/2 (for g — 1 these are n x w 0.3771 and 
n x ps 0.6229). Inside the intermediate phase these two 
peaks would move towards each other (and n x = 1 /2) as 
the dot is made more isolated, eventually coalescing as 
v approaches v C 2- An alternative way of observing the 
existence of the half-integer (and also integer) plateaus 
is through detecting the associated hysteresis phenomena 
which must accompany the first-order population jumps. 
We believe that the observation of hysteresis phenom- 
ena is a good gauge to distinguish between the predicted 



continuous and discontinuous n(n x ) dependencies. 

We believe, additionally, that these plateaus should 
be accessible numerically. Monte Carlo simulations of 
actions related to that found in (|]) have already been 
successively accomplished. Specifically setting b = in 
(||) gives us an action appropriate to describing tunneling 
between quantum Hall edges. Monte Carlo simulations 
of this action were done in pdj| which later were shown to 
match exact scaling curves [fl2[. Moreover very accurate 
Monte Carlo simulation of the classical one-dimensional 
Ising model with interactions decaying as inverse-square 
of distance was recently completed |50| , and the results 
were found to be in very good agreement with theoreti- 
cal predictions |2§| , ^9[| . We are thus currently examining 
numerically the possibility of observing the half- integer 
plateaus |43|| . 

Having outlined the prospects of observing the inter- 
mediate phase, we may note that in the context of ad- 
sorption phenomena the existence of a regime with stable 
fractional film thickness centered about half-layer posi- 
tion is experimentally well-established H|. In fact, the 
effects observed are even more puzzling - in a range of 
parameters variations in the external pressure lead to 
film thickness jumps only between half-integer layer cov- 
erages. The regime analogous to our finding - integer- 
half-integer adsorption curve was discussed theoretically 
p4j , but not yet seen experimentally. 

We would also like to mention a recent Little-Parks 
type experiment on small disordered AuZn cylindrical 
films where a resistance oscillation with a period of 
$o/2, half of the flux quantum, was observed p5| . The 
new resistance minima develop below certain tempera- 
tures and coexist with the standard Little-Parks min- 
ima corresponding to the integer number of flux quanta. 
This closely resembles our prediction for the intermediate 
phase translated into the superconductivity language. 

The persistence of the half-integer minimum in the 
population dependence of the ground-state energy at the 
charge degeneracy points is the ultimate reason why the 
physics of the Kondo/Ising systems is different from that 
of the tricritical systems. The presence of an extra min- 
imum in the latter case makes it more difficult for the 
electrons to enter the dot, thus enhancing the effects of 
Coulomb blockade. 

Putting aside the existence of the intermediate phase, 
in the g < 1 regime there is nonetheless a qualitative 
correspondence between the results of Sections III and 
IV - for fixed g and sufficiently small v (sufficiently large 
bare tunneling amplitude Ao) the n(n x ) dependence is 
continuous while above a critical v (below critical Ao) 
the n(n x ) dependence is discontinuous - it is a modified 
staircase with plateaus centered around integer dot pop- 
ulation. The correspondence can be made more quan- 
titative by comparing the dot population jump between 
two neighboring integer plateaus numerically evaluated 
at the upper critical curve v c i from (0) - (O) with the 
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universal prediction An = g 1 / 2 for the two-state sys- 
tem (i.e. a ferromagnetic Ising model with interactions 
decaying as inverse square of distance). We found that 
whenever the width of the intermediate phase is small, 
i.e. for (i) < g < 1 and b — ► or (ii) g — * and arbi- 
trary b, the population jump computed from ( pd] ) - fll3| ) 
is well approximated by g 1 ^ 2 . For g — ► 0, a little algebra 
shows that An ~ g 1 / 2 is a consequence of Eqs. ( |l2| ) and 
(|l3|), and the prefactor can be established numerically. 
We find that An = 1.019913 1 / 2 as g — ► which is very 
close to the theoretical expectation An = g 1 ^ 2 . 

For the experimentally most significant case when the 
electrons in the reservoirs are noninteracting, g = 1, the 
discrepancy between the systems of the Kondo/Ising and 
the tricritical subclasses is largest. In the former case a 
continuous n(n x ) dependence is predicted for arbitrarily 
small bare tunneling amplitude, while in the latter 
continuous, intermediate staircase, and modified integer 
staircase are predicted, as in Fig. 2, as the tunneling am- 
plitude decreases (v increases). This difference between 
the two types of behavior is unsurprising as the treatment 
of Section III for the case g = 1 was already extremely 
delicate - the slightest change in the physics such as a 
stronger suppression of the number fluctuations on the 
dot might change the outcome qualitatively. Experimen- 
tal or numerical observation for g = 1 of the modified 
Coulomb staircase with first-order jumps between inte- 
ger population plateaus will provide indirect evidence for 
the intermediate phase - they either exist together (tri- 
critical type) or do not exist at all (Kondo/Ising type). 
It would be also interesting to verify the universality of 
the n(n x ) dependence ( fl^ ) by varying the size and shape 
of large dots. 

Beyond the existence of the half-integer plateaus, an 
important issue arises in relating the smearing of the 
Coulomb staircase found in M| with our results. One 
pos sible answer is that the effective action implicit in 
fLOf belongs to what we call the Kondo/Ising subclass 
(specifically a Kondo Hamiltonian with the coupling J z 
set to 0) while our variational finding of a staircase oc- 
curs solely for actions falling in the tricritical subclass. 
In Section III we demonstrated that after the dot prob- 
lem is truncated to a two-state system, any model with 
an effective free-fermion reservoir will predict a smeared 
Coulomb staircase. Because in the classical limit of a 
closed dot the free-fermion case is marginal, it is also 
important to critically reexamine the validity of the two- 
state approximation. Let us assume that the gate volt- 
age is tuned at a charge degeneracy point, and that the 
full bare potential (|^) has metastable minima ignored in 
the two-state approximation. Because these extra min- 
ima can now be populated with small finite probability, 
the full potential (^) suppresses electron number fluctu- 
ations on the dot somewhat stronger than its two-state 
truncation discussed in Section III [^6). Whether this 
suppression is strong enough to restore a true staircase, 



remains an open problem. 

We have also conducted a study of the spinful problem 
akin to the variational analysis of Section IV. The results 
will be presented elsewhere. Similar to the spinless case, 
the spin- 1/2 Coulomb blockade problem is dual to that 
of a two-junction SQUID pjfl . Although this problem is 
more involved, for the case when the electrons are non- 
interacting, the model is a natural generalization of (Q). 
We find that for sufficiently small transparency of the 
constriction there are two staircase phases in close corre- 
spondence with the spinless case. Experimentally p"5[ for 
spinful electrons, the dot population was found to change 
continuously with the gate voltage for a nearly open con- 
tact while developing smeared step-like structures with 
the contact nearly closed. The shape of the smeared step 
for the almost closed dot was found to be consistent with 
a finite temperature generalization of the perturbation 
theory done in |10|. However the temperature of the 
dot was far above the Kondo temperature found in [fl0|| , 
the scale governing the putative smearing of the steps at 
zero temperature. Whether the steps remain smeared at 
T = is thus an unanswered experimental question. 

We hope that experimental attempts will be made in 
the near future to test our predictions. 
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